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ABSTRACT 

I present ignition models for Type I X-ray bursts and superbursts from the ultracompact binary 4U 1820-30. 
A pure helium secondary is usually assumed for this system (which has an orbital period re 11 min), however 
some evolutionary models predict a small amount of hydrogen in the accreted material (mass fraction X ~ 0.1). I 
show that the presence of hydrogen significantly affects the Type I burst recurrence time if X > 0.03 and the CNO 
mass fraction >3x 10~ 3 . When regularly bursting, the predicted burst properties agree well with observations. 
The observed 2—4 hour recurrence times are reproduced for a pure He companion if the time-averaged accretion 
rate is (M) re 7-10 x 10" 9 M Q yr" 1 , or a hydrogen-poor companion if (M) re 4-6 x 10" 9 M Q yr" 1 . This result 
provides a new constraint on evolutionary models. The burst energetics are consistent with complete burning, and 
spreading of the accreted fuel over the whole stellar surface. Models with hydrogen predict ~ 10% variations in 
burst fluence with recurrence time, which perhaps could distinguish the different evolutionary scenarios. I use 
the M's determined by matching the Type I burst recurrence times to predict superburst properties. The expected 
recurrence times are re 1-2 years for pure He accretion (much less than found by Strohmayer and Brown), and 
re 5-10 years if hydrogen is present. Determination of the superburst recurrence time would strongly constrain 
the local accretion rate and thermal structure of the neutron star. 

Subject headings: accretion, accretion disks — X-rays:bursts — stars:neutron — stars individual (4U 1820-30) 



1. INTRODUCTION 

Type I X-ray bursts are well-understood as thermonuclear 
flashes on the surface of an accreting neutron star (see Lewin, 
van Paradijs, & Taam 1995; Strohmayer & Bildsten 2003 for re- 
views). The basic theory was outlined more than twenty years 
ago (for an overview, see Fujimoto, Hanawa, & Miyaji 1981; 
Bildsten 1998), successfully explaining the typical burst ener- 
gies (10 39 -10 40 ergs), durations (~ 10-100 s), and recurrence 
times (hours to days). However, detailed comparisons of theory 
and observations have had mixed success (Fujimoto et al. 1987; 
van Paradijs, Penninx, & Lewin 1988; Bildsten 2000; Galloway 
et al. 2003). 

The ultracompact binary 4U 1820-30 (P orb = 1 1 .4 min; Stella, 
Priedhorsky, & White 1987) is a promising system for such 
a comparison. It has a known distance, being located in the 
metal-rich globular cluster NGC 6624 ([Fe/H] re -0.4, dis- 
tance 7.6 ±0.4 kpc); Rich, Minniti, & Liebert 1993; Kuulk- 
ers et al. 2003). It undergoes a regular re 176 day accre- 
tion cycle (Priedhorsky & Terrell 1984), switching between 
high and low states differing by a factor re 3 in accretion rate 
(e.g. see the RXTEI 'ASM light curve in Figure 1 of Strohmayer 
& Brown 2002). In the low state, regular Type I bursts are seen 
(±23 days around the minimum luminosity; Chou & Grindlay 
2001), with recurrence times re 2-4 hours (Clark et al. 1976; 
Clark et al. 1977; Haberl et al. 1987; Cornelisse et al. 2003). 
No bursts are seen during the rest of the cycle, implying a 
"non-bursting" mode of burning (e.g. Bildsten 1995). In the 
low state, 4U 1820-30 has also shown an extremely energetic 
(~ 10 42 erg) "superburst", likely due to deep ignition of a car- 
bon layer (Strohmayer & Brown 2002; see Kuulkers et al. 2002 
and Strohmayer & Bildsten 2003 for a summary of superburst 
properties). 

Different evolutionary scenarios for 4U 1820-30 have been 
proposed. Those involving direct collision of a red giant and a 
neutron star (Verbunt 1987), or formation of a neutron star-main 
sequence binary by tidal capture (Bailyn & Grindlay 1987) or 



an exchange interaction (Rasio, Pfahl, & Rappaport 2000) fol- 
lowed by a common envelope phase, lead to accretion of pure 
helium from the helium core of the red giant. A different picture 
is that the mass transfer starts just after central hydrogen ex- 
haustion (Tutukov et al. 1987; Fedorova & Ergma 1989; Podsi- 
adlowski, Rappaport, & Pfahl 2002), in which case the accreted 
material is again mostly helium, but contains some hydrogen 
(re 5-35% by mass; Fedorova & Ergma 1989; Podsiadlowski, 
Rappaport, & Pfahl 2002). 

The observed energetics and recurrence times of the Type 
I bursts from 4U 1820-30 fit well with a picture of accumu- 
lation and burning of helium-rich material. The ratio of per- 
sistent fluence (persistent flux integrated over the burst recur- 
rence time) to burst fluence is a re 120 (Haberl et al. 1987). 
For a gravitational energy GM/R re 200 MeV per nucleon, 
this implies a nuclear energy release re 1.6 MeV per nu- 
cleon, as expected for helium burning to iron group elements. 
The mass of helium required to power the observed energet- 
ics is roughly the mass accreted in the burst recurrence time 
at the inferred accretion rate. The burst fluence re 3.5 x 
10~ 7 erg cm" 2 (Haberl et al. 1987) implies a total burst energy 
2.5 x 10 39 erg (d/7.6 kpc) 2 , or a mass of helium AM re 1.6 x 
10 21 g. The persistent luminosity when bursts are seen is Lx re 
2.8 x 10 37 erg s" 1 {d/7.6 kpc) 2 (F x /4 x 10" 9 erg cm" 2 s" 1 ), giv- 
ing an accretion rate M re 1 .5 x 10 17 g s" 1 re 2.4 x 10" 9 M© yr" 1 
for a 1.4 Mq neutron star with radius R = 10 km. At this rate, 
a mass AM is accreted in 3 hours, in excellent agreement with 
the observed recurrence times. 

A previous comparison of theory with observations of 
4U 1820-30 was carried out by Bildsten (1995, hereafter B95), 
who conducted time-dependent simulations of pure helium 
burning on accreting neutron stars. He found good agreement 
with observed energetics and recurrence times, but for a some- 
what hotter base temperature than expected. In this paper, I 
calculate ignition models for Type I bursts from 4U 1820-30, 
and survey the conditions necessary to ignite helium at the in- 
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ferred mass AM. I allow for a small amount of hydrogen in the 
accreted material, and show that this has a large effect on the 
burst recurrence time, because of extra heating from hot CNO 
hydrogen burning as the layer accumulates. Finally, I discuss 
what we might learn from simultaneous modelling of Type I 
bursts and superbursts. The plan of the paper is as follows. The 
main results are presented in §2. I calculate the critical mass 
needed for ignition, and show the effect of hydrogen on recur- 
rence times and energetics. I compare these results to observa- 
tions of Type I bursts, and to evolutionary models. Finally, I use 
the constraints from Type I bursts to construct ignition models 
for superbursts from 4U 1820-30, and compare to the model 
of SB02. I summarize in §3, and discuss the remaining open 
issues. 

2. IGNITION MODELS AND COMPARISON TO OBSERVATIONS 

2.1. Calculation of Ignition Conditions 

The physics and different regimes of nuclear burning on ac- 
creting neutron stars have been described by several authors 
(for reviews, see Lewin, van Paradijs, & Taam 1995; Bildsten 
1998). Here, I follow the X-ray burst ignition calculations of 
Cumming & Bildsten (2000, hereafter CB00). I adopt plane- 
parallel coordinates in the thin accreted layer, and work in terms 
of the column depth y = P/g, where dy = -pdz, P is the hydro- 
static pressure, p is the density, and g the surface gravity. I 
assume a 1.4 M Q neutron star with radius R = 10 km, giving 
g = (1 +z)GM/R 2 = 2.45 x 10 14 cm s" 2 where the gravitational 
redshiftisz = l-(l-2GM//?c 2 r 1 / 2 = 0.31 (similar to the value 
« 0.35 recently inferred from absorption lines in XMM burst 
spectra of EXO 0748-676; Cottam, Paerels, & Mendez 2002). I 
include general relativistic effects when calculating observable 
quantities. Cumming et al. (2002) showed that the Newtonian 
equations of CB00 have the same form under general relativity, 
as long as y and m refer to rest mass rather than gravitational 
mass. 

The thermal profile of the accumulating layer is described by 
the heat equation F = (4acT 3 /3n)(dT /dy), where the opacity k 
is calculated as described by Schatz et al. (1999), and the en- 
tropy equation dF/dy = -e, where F is the heat flux, and e con- 
tains contributions from hot CNO hydrogen burning and com- 
pressional heating. Compressional heating always plays a mi- 
nor role at these accretion rates, giving <~ cpT « 5ksT /2fim p « 
0.02 7g MeV per nucleon. The hot CNO energy production rate 
is e H = 5.8 x 10 13 erg g" 1 s" 1 (Z/0.01) (Hoyle & Fowler 1965), 
where Z is the CNO mass fraction. I follow the changing hydro- 
gen abundance with depth as described in CB00 (see eqn. [3] of 
that paper). Unless otherwise stated, I set Z = 0.01. Triple al- 
pha burning during accumulation is not included, this has only a 
small effect on the ignition conditions. I assume the individual 
CNO abundances are those of the equilibrium hot CNO cycle 
( 14 and 15 0). In the pure He case, the accreted 14 N survives 
to a density w 10 6 g cm" 3 , when it converts to 14 C by electron 
capture (Hashimoto et al. 1986). For the conditions here, no 
further processing occurs until helium ignites via triple alpha. 

To find the thermal profile, I repeatedly integrate downwards 
through the accreted layer, iterating until F = F cmst at the base, 
where F crust = Q cmst (m) is the heat flux from deeper layers. 
Here, Qcmst is the energy per nucleon released in the crust 
from pycnonuclear and electron capture reactions which es- 
capes from the surface. Brown (2000) showed that Q C rust 

«0.1 

MeV per nucleon « 10 17 erg g" 1 for m > 10 4 g cm" 2 s" 1 . I as- 
sume F cmst is set by the time-averaged local accretion rate (m), 



rather than the instantaneous rate m, because the thermal time in 
the crust is much longer than the w 6 month luminosity cycle. 
Throughout this paper, I assume the instantaneous rate when 
bursts are seen is half the time-average rate, m = (rh) /2 (as in- 
ferred from the RXTE/ ASM lightcurve; e.g. Figure 1 of SB02). 

The ignition column depth y\ gn is found by increasing the 
depth of the layer until the condition for ignition, 

"^3a _ ^^cool n. 

~df~~ dT 

(Fujimoto, Hanawa, & Miyaji 1981; Fushiki & Lamb 1987b; 
Bildsten 1998), is met at the base. Here is the triple 
alpha energy production rate (Fushiki & Lamb 1987a), and 
e cooi = acT 4 /3ny is a local approximation to the radiative cool- 
ing. The ignition criterion (1) says that ignition occurs when the 
heating from triple alpha is more temperature sensitive than the 
cooling, and is large enough to change the temperature of the 
layer. Following CB00, we include a multiplicative factor of 
1 .9 in £3 Q to account for proton captures on 12 C when hydrogen 
is present. 

The advantage of this "semi-analytic" approach is that it al- 
lows a survey of parameter space, while still giving a good esti- 
mate of the ignition depth. CB00 compared the ignition condi- 
tions predicted by equation (1) with time-dependent models in 
the literature. For mixed H/He ignition, the agreement was to 
10-30%, and better than a factor of two for pure He ignitions. 
In addition, I have compared results for pure helium with those 
of B95. For a given m, I adjust the flux from the base to match 
the base temperature specified by B95. For the m's relevant to 
this paper, I find that the ignition column depth from equation 
(1) underpredicts the ignition column by 10-30%. 

I estimate the burst energy by assuming complete burning 
of the accreted layer, and taking a nuclear energy release Q mc = 
1 .6+4.0 (X) MeV per nucleon, where (X) is the mass-weighted 
mean X in the layer at ignition. This expression for Q nuc in- 
cludes 35% energy loss in neutrinos during hydrogen burning 
(Fujimoto et al. 1987). This is appropriate for rp-process hydro- 
gen burning up to iron group, and may be an overestimate for 
small Xq, in which case we underestimate the burst energy by up 
to « 10% for (X) ^0.1. Assuming the layer of fuel covers the 
whole star, the total burst energy is Eburst = 4nR 2 yi sn Q nuc / ( 1 +z), 
where the (1 + z) factor accounts for gravitational redshift. 

2.2. Effect of Hydrogen 

Hydrogen burning gives w 7 MeV per nucleon, so that 
even a small amount of hydrogen burning during accumu- 
lation dominates the « 0.1 MeV per nucleon coming from 
the crust. For pure helium accretion, f crust = (m) Qcmst « 
10 21 erg cm 2 s" 1 (m} 4 gn sets the temperature of the accumu- 
lating layer. If hydrogen is present, the flux from the hot CNO 
cycle is F H = enyH = 5.8 x 10 21 erg cm 2 s" 1 yn,$ (Z/0.01). Here, 
yu is the column density of the layer which is burning hydro- 
gen: this is either the column depth at which hydrogen runs 
out, or the depth at which helium ignites, whichever occurs 
first. Whether F H > F crust depends on y H since this measures 
the amount of hydrogen burning. The time to burn all the hy- 
drogen is t H = E H X / e H , or 

-"-(^(fr). 

where E H = 6.0 x 10 18 erg g 1 is the energy release per gram 
from hydrogen burning (unlike CB00, we include the energy 
release from neutrinos here as given by Wallace & Woosley 
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1981), and Xq is the initial hydrogen abundance. If the burst 
recurrence time f recur < tn, then helium ignition occurs be- 
fore hydrogen exhaustion, giving y H = rht recul . Then requiring 
Fh > -Fcrust gives 

Z>3.5xl0- 3 e, 7 (^). (3) 

If frecur > fe, all of the hydrogen is burned before helium ignites, 
giving y H = mt H = 1.2 x 10 8 g cm" 2 (X /0.1)(m 4 /1.2)(0.01/Z). 
Then 

X > 0.034 Qn(J£r) (4) 

is required for F H > F crus t- 

Equations (3) and (4) show that if hydrogen burning is to 
affect the temperature profile in the layer, both the amount of 
hydrogen and the metallicity must be large enough. Evolu- 
tionary models for 4U 1820-30 with hydrogen give X > 0.1 
(e.g. Podsiadlowski et al. 2002), satisfying equation (4). In ad- 
dition, NGC 6624 is metal-rich ([Fe/H] « -0.4; Rich, Minniti, 
& Liebert 1993) so that the accreted metallicity satisfies equa- 
tion (3) (although spallation in the accretion shock could re- 
duce the metallicity in the accreted layer; Bildsten, Salpeter, & 
Wasserman 1992). 

2.3. Comparison with Observations 

Several authors have presented observations of regular burst- 
ing from 4U 1820-30 (the first source to exhibit X-ray bursts; 
Grindlay et al. 1976). Clark et al. (1976) saw 10 bursts with 
SAS-3 with a recurrence time of 4.4 ± 0. 17 h (a 4% variation), 
and mean fluence « 2.9 x 10~ 7 erg cm" 2 (1-10 keV). The per- 
sistent flux was « 1/5 of the peak flux seen in the high state. 
Further SAS-3 observations of the transition from low to high 
state found 22 bursts (Clark et al. 1977). The recurrence time 
decreased from 3.4-2.2 hours as the flux increased by a factor 
of w 5. The mean burst fluence was « 4 x 10" 7 erg cm" 2 . The 
decay time of the bursts decreased slightly as flux increased. 
Haberl et al. (1987) saw 7 bursts in a 20 hour EXOSAT obser- 
vation with a recurrence time of 3.2 ± 0.05 hours, mean fluence 
(3.48 ± 0. 16) x 10" 7 erg cm" 2 , and a = 142 ± 10. The persis- 
tent flux was « 4.3 x 10~ 9 erg cm" 2 s" 1 . Most recently, Bep- 
poSAX/WFC accumulated « 7 Ms of exposure of the Galac- 
tic Center region, observing 49 bursts from 4U 1820-30 (Cor- 
nelisse et al. 2003; Kuulkers et al. 2003). These bursts show 
a transition from regular to irregular bursting as persistent flux 
increases (Cornelisse et al. 2003). 

The simplest model is to consider pure helium accretion 
at the rate inferred from the X-ray luminosity in §1, and to 
take gcrust = 0.1 MeV per nucleon, as found by Brown (2000). 
As a useful reference point, I refer to this accretion rate as 
nix = 1 .2 x 10 4 g cm" 2 s" 1 . Figure 1 shows the temperature pro- 
file at ignition for this case (lower curve). The recurrence time 
and burst energy are w 5 times larger than the observed values. 
Ignition occurs at a column depth « 7 x 10 8 g cm" 2 , giving a 
recurrence time of 16.5 h, and burst energy 1.1 x 10 40 erg. A 
model with hydrogen gives much better agreement. The up- 
per curve in figure 1 shows the same model, but with Xq = 0.1. 
The extra heating from hot CNO burning leads to earlier ig- 
nition, y w 2 x 10 8 g cm" 2 , giving f recur = 4.4 h and £t,urst = 
3.1xl0 39 erg. 

Unfortunately, this result can not be taken as evidence for the 
presence of hydrogen. A larger flux from below will also heat 
the layer and lead to earlier ignition. For pure helium accretion 



at hi = m x , I find that Q cmsl = 0.4 MeV per nucleon is required 
to bring the recurrence time down to w 3 hours. This is larger 
than expected from heating in the crust (Brown 2000). How- 
ever, additional heating might be provided by a previous X-ray 
burst. Most of the energy released in a burst is radiated from the 
surface of the star; however, some energy is conducted inwards, 
and emerges on a longer timescale. This effect has been empha- 
sized by Taam and coworkers (Taam 1980; Taam et al. 1993), 
who refer to it as "thermal inertia". To estimate the size of 
this effect, I assume a fraction / of the nuclear energy flows in- 
wards, and is released at a constant rate between bursts. The ex- 
tra flux is then fE\, urst /4-TrR 2 t recm , giving a contribution to Q CIUSt 
of w 0.1 MeV per nucleon (//0.1)(3 h/t Kcm )(E 39 /3)(m x /m). 
This is likely an overestimate, since the energy release de- 
creases over time. A better estimate requires time-dependent 
simulations of a sequence of bursts. For now, I incorporate this 
uncertainty into g crust . 

An additional uncertainty is the Lx-rh relation, which allows 
some freedom in setting rh. Pure helium models nicely agree 
with observations if hi w 1.5-2 ni x (for Q cmst = 0.1-0.2 MeV 
per nucleon). This is shown in Figure 2, in which f recur and 
£burst are plotted against ni. For each choice of ni, the time- 
averaged rate is assumed to be (rh) = 2m. The vertical bars at 
m = m x show the observed range of recurrence time and energy. 
As we found above, hydrogen models require hi w ni x to give 
agreement. 

Table 1 shows the same result for some specific models. 
Here, I choose Xq, Z, and Qcmst, and then vary m to match 
the recurrence time seen by Haberl et al. (1987), f recur = 3.2 
hours. There a few points worth noting. First, the ignition 
depth in the helium models is sensitive to (2 crust , since this sets 
the temperature of the accumulating layer. As Q cmst increases, 
y; gn decreases. Second, hydrogen models show two regimes of 
behavior. If f recur < tu (see eq. [2]), hydrogen burning occurs 
throughout the accumulating layer. Then y; gn depends on the 
metallicity Z which determines the amount of hot CNO burn- 
ing, but not on gcrust or rh. As Z decreases, y\ gn increases. An 
extreme case is Model 7 in Table 1 (Z = 10" 3 ), in which case lit- 
tle heating occurs and the ignition conditions are similar to pure 
helium. If f recur > tu, hydrogen runs out and ignition occurs in 
a pure helium layer. Then y\ gn is sensitive to rh, explaining the 
steep increases in f recur and £burst for the H models in Figure 2 
(see CB00 for more discussion of these different regimes). 

In general, the burst energy is overpredicted in these mod- 
els, by factors of up to « 70%. This discrepancy is easily ac- 
counted for by anisotropic burst emission (e.g. Lapidus & Sun- 
yaev 1985; Fujimoto 1988) (although there may be other ex- 
planations, see discussion in §3.1). For the remainder of this 
section, I renormalize the predicted fluence to agree with the 
EXOSAT mean value (Haberl et al. 1987). 

Figure 3 shows the expected variations in burst fluence and 
recurrence time as a function of persistent flux. I show hydro- 
gen models with metallicity Z = 0.01 (X = 0.1 and X = 0.2), a 
low metallicity model Z = 10~ 3 (X = 0.1), and a helium model 
with g crust = 0.1. I assume g crust is independent of accre- 
tion rate, not including any variations with hi because of ther- 
mal inertia. In each case, I normalize the F x -ni relation to 
give f recur = 3.2 h for F x = 4.3 x 10~ 9 erg cm" 2 s" 1 (Haberl et 
al. 1987). To do this, I write F x = L x /4nd 2 = mQ gtm (R / d)% 1 , 
where g grav = c 2 z/(l+z) is the gravitational energy release (221 
MeV per nucleon for z = 0.31), and the factor £ p accounts for 
the anisotropy of the X-ray emission, as well as additional un- 
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certainties such as luminosity emerging at other wavelengths. 
Given the m for the model with ? recur = 3.2 h, I set 

r i( R/d V_ 1-1 (Fx,-9\ / 221MeV \ 
^ V10km@7.6kpcy m 4 \ 4.3 /V Ggrav / 

which gives the correct normalization. The values of 
^(Re/dj.e) 2 implied are 0.5 (0.8) for the pure He (X = 0.1) 
model. Since the thermal time in the crust is much longer than 
the 6 month luminosity cycle, I keep (m) fixed as F x varies, 
at the value appropriate for f recur = 3.2 h. 

The expected variation of recurrence time with m agrees 
qualitatively with the observed behavior in that recurrence time 
drops with increasing Fx. Clark et al. (1976) found f recur = 
4.4 when the flux was « 1/5 of its peak value. Estimating 
the peak flux to be roughly 3 times the EXOSAT value, or 
« 1.2 x 10~ 8 erg cm" 2 s" 1 , gives Fx = 2.4 x 10~ 9 erg cm" 2 s" 1 
for these observations, less than the EXOSAT value and con- 
sistent with a longer f recU r- However, the dependence of f recur on 
m observed by Clark et al. (1977) is not as steep as oc 1/m pre- 
dicted in the models. They report that f recur decreased from 3.4 
hours to 2.2 hours as the persistent flux increased by a factor 
of 5. In the models, only the X = 0.1 model shows a signif- 
icant deviation from a 1/m law, but in the opposite direction. 
For f recur < t H = 3 hours, f recur falls off more steeply than 1/m, 
opposite to the observed trend. 

In Figure 4, I plot fluence against / recur directly. The flu- 
ence is renormalized to agree with the mean EXOSAT value, 
highlighting the variation in fluence with f recur . Hydrogen mod- 
els show a change in fluence of 10% as the recurrence time 
varies over the observed range, whereas the helium and low 
metallicity models have roughly constant fluence. When hy- 
drogen is present, the variation arises because the composition 
of the layer at ignition changes according to the amount of hy- 
drogen burning that occurs during accumulation. I also plot the 
observed values in Figure 4, showing the range of recurrence 
times for the Clark et al. (1977) observations, measurements 
for individual bursts from Haberl et al. (1987), and the mean 
fluence reported by Clark et al. (1976). There is a m 30% vari- 
ation in the reported fluences. However, the error in the SAS-3 
fluence measurements is likely larger than this. For example, 
a reanalysis by Clark (reported by Vacca et al. 1986) gave a 
« 40% reduction in the peak luminosities of individual bursts. 
Therefore any conclusions drawn from Figure 4 are probably 
unwarranted. The comparison is complicated by intrinsic scat- 
ter in the burst fluences. The scatter in the Haberl et al. (1987) 
points in Figure 4 is 4%, comparable to the predicted variation 
with recurrence time if hydrogen is present. Nonetheless, look- 
ing for variations in fluence with recurrence time or persistent 
flux is a promising way to distinguish between models with and 
without hydrogen. It perhaps could be addressed with the Bep- 
poSAX data (Cornelisse et al. 2003). 

2.4. Comparison with Evolutionary Models 

The requirement that m « m x if hydrogen is present, or 
m « \.5-2m x for pure helium, is a new constraint on evo- 
lutionary models. Evolutionary models for 4U 1820-30 pre- 
dict values of M and Xq (Fedorova & Ergma 1989; Podsiad- 
lowski et al. 2002). Table 2 gives the resulting ignition condi- 
tions for Type I bursts. There is a remarkable agreement be- 
tween the m's required to match Type I burst properties, and 
those predicted by evolutionary models (although it should be 
noted that these predicted m's are rather uncertain, e.g. see 



Podsiadlowski et al. 2002). For example, if the companion 
is a cold helium white dwarf, then gravitational radiation an- 
gular momentum losses give M CR « 8.8 x 10~ 9 M Q yr" 1 (for 
a companion mass of 0.068 M ; Podsiadlowski et al. 2000), 
or m :=y (m) /2 w 2.2 x 10 4 g cm" 2 s" 1 (R/10 km)" 2 , consistent 
with the values we find in Table 1 for X = 0. In addition, evo- 
lutionary models with hydrogen tend to have lower accretion 
rates than pure helium models, by about a factor of two, so that 
in general both are consistent with the data. There are partic- 
ular models, however, that do not give good agreement. For 
example, Podsiadlowski et al. (2002) give a hydrogen model 
with M = 1.1 x 10~ 8 M C3 yr" 1 , which has a recurrence time of 
1 .4 hours, outside the observed range. 

2.5. Superburst Ignition Conditions 

A superburst was seen from 4U 1820-30 on 1999 Sep 9 
(Strohmayer & Brown 2002). The source was in the low state, 
with F x = 3.5 x 10" 9 erg cm" 2 s" 1 . The total fluence of the su- 
perburst was 2.7 x 10" 4 erg cm" 2 , giving a total energy release 
of « 2 x 10 42 ergs (almost a factor of 1000 larger than a normal 
Type I burst). SB02 modelled the superburst from 4U 1820- 
30, proposing that carbon production during stable burning in 
the high state leads to a thick carbon/iron layer (mass w 10 26 g, 
carbon mass fraction X c > 0.3) which ignites to give the super- 
burst. 

SB02 estimated the accretion rate for their model from the X- 
ray luminosity, finding m w m X - I now calculate ignition con- 
ditions for the 4U 1820-30 superburst using the m's inferred 
from Type I bursts. I follow the carbon ignition calculations of 
Cumming and Bildsten (2001, herafter CB01), except I take the 
heavy element to be 56 Fe (CB01 considered much heavier ele- 
ments than iron, made during rp process burning of hydrogen). 
I integrate downwards through the carbon/iron layer, assuming 
a constant heat flux from the crust F crust = Q cm st (m) (I assume 
that the heat flux is the same as that heating the helium layer 
between normal Type I bursts). I calculate the carbon ignition 
depth using a similar criterion to equation (1). To estimate the 
superburst energy, I assume that carbon burning gives 1 MeV 
per nucleon. The amount of carbon is uncertain, and depends 
on the fraction of stable H/He burning, and the carbon yield. To 
allow a direct comparison with SB02, 1 assume that the carbon 
mass fraction Xc is 30% in these models. 

Table 3 gives the resulting ignition conditions. The model 
numbers correspond to those in Table 1, except for the model 
labelled "SB02", which has m = m X - This gives f recur = 14 years, 
and total energy release w4 x 10 43 ergs, similar to the values 
found by SB02. For the helium models, the larger m than as- 
sumed by SB02 gives an order of magnitude shorter recurrence 
time, w 1-2 years rather than 14 years (this sensitivity of re- 
currence time to m was noted by SB02, see also Figure 3 of 
CB01). For the X = 0.1 models, the recurrence time is longer, 
4-10 years, depending on gcrust- This shows that the superburst 
recurrence time, currently not well constrained, is a powerful 
discriminant between different models. 

Note that the total nuclear energy release can be much greater 
than the observed superburst energy if neutrino emission plays 
an important role. This was emphasized by SB02 for their 
model of the 4U 1820-30 superburst, in which only w 10 42 erg 
(from a total of 4 x 10 43 ergs) is radiated from the surface. How- 
ever, the amount of neutrino emission is strongly dependent on 
the peak temperature, and therefore carbon fraction and ignition 
depth (CB01). Time-dependent models are therefore necessary 
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to predict the superburst fluence. 

CB01 showed that for small carbon mass fractions, the car- 
bon burns stably as the layer accumulates. Table 3 gives the 
minimum Xq required for unstable ignition for each model. The 
range of values is Xc, m m = 0.08-0.2. Less carbon is required 
for smaller g crust , as shown by CB01. The values of Xc,mm 
given here are a little larger than those shown in CBOl's Fig- 
ure 2 because the less massive heavy element ( 56 Fe rather than 
104 Ru) results in a shallower temperature gradient. The recur- 
rence time for Xq = -Xc.min is « 10% lower than the recurrence 
time forX c = 0.3. 

3. SUMMARY AND DISCUSSION 

3.1. Type I bursts 

In general, the models of Type I X-ray bursts presented in §2 
give good agreement with observed recurrence times and en- 
ergetics. The observed 2-4h recurrence times are reproduced 
for a pure He companion if the time-averaged accretion rate is 
(M) w 7-10 x 10~ 9 Mq yr _1 (R/\0 km) 2 , or a hydrogen-poor 
companion if (M) «4-6 x 10- 9 M Q yr" 1 (fl/10km) 2 . This dif- 
ference comes about because of extra heating from hot CNO 
burning during accumulation when hydrogen is present, lead- 
ing to earlier ignition, and requiring a lower m to match the 
observed recurrence times. For hydrogen to affect the temper- 
ature profile in the layer requires X > 0.03 and Z > 3 x 10" 3 . 
If spallation reduces the metallicity in the accreted layer (Bild- 
sten, Salpeter, & Wasserman 1992), hydrogen burning will no 
longer play a role, in which case the accreting layer behaves 
like pure helium. 

These conclusions constrain evolutionary models, although 
both pure helium evolutionary models and those with hydrogen 
can be found that are consistent with the data. A promising way 
to distinguish between different evolutionary scenarios is accu- 
rate fluence measurements as a function of m. Models with hy- 
drogen predict « 10% variations in burst fluence as recurrence 
time varies, whereas pure helium models show a much smaller 
variation. 

For most models, recurrence time is expected to vary as 
w l/m, because the ignition depth (determined by hot CNO 
burning or by the flux from the crust) is insensitive to accre- 
tion rate. However, SAS-3 observations showed a much shal- 
lower dependence than the expected l/m: Clark et al. (1977) 
observed the persistent flux increase by a factor of « 5 while 
frecur went from 3.4 to 2.2 hours. Studies of timing and spectral 
behavior of LMXBs indicate that Lx is not necessarily well- 
correlated with the underlying accretion rate (e.g. van der Klis 
et al. 1990). The observed behavior then requires m to increase 
less quickly than L x . It would be interesting to study burst prop- 
erties as a function of an accretion rate indicator (e.g. position 
in the color-color diagram) rather than Lx . Variations in heating 
with m due to thermal inertia might play some role, something 
that should be checked with time-dependent simulations. Al- 
ternatively, these observations might indicate new physics, for 
example that the accreting material covers only part of the star 
at lower accretion rates, reducing the amount by which the lo- 
cal accretion rate changes (a possibility suggested for hydrogen 
accretors by Bildsten 2000). An interesting exception to the 
predicted l/m behavior is for models with Xq « 0.1, when the 
time to burn all the hydrogen is « 3 hours, in the range of ob- 
served recurrence times. ForXo ~ 0.1 and f recur > 3 hours, the 
recurrence time is predicted to increase rapidly with decreasing 
m, because hydrogen burning no longer occurs throughout the 



layer. However, this is even less consistent with the observed 
variations. 

The burst energy is overpredicted in these models, in com- 
mon with studies of other bursters. Galloway et al. (2003) find 
a factor of 3-4 discrepancy for GS 1 826-24, and following Bild- 
sten (2000) propose that the accreted fuel only partially covers 
the neutron star surface. For4U 1636-56, Fujimoto et al. (1987) 
proposed that vertical mixing driven by rotational shear in the 
incoming material led to early ignition, and reduced burst en- 
ergies. Here, I find discrepancies of less than a factor of 2 
for 4U 1820-30. This is easily accomodated, for example, by 
anisotropic burst emission (Lapidus & Sunyaev 1985; Fujimoto 
1988), and does not require incomplete spreading of the fuel or 
extra mixing. Note however that the discrepancy increases if 
the neutron star radius is larger than the 10 km assumed here 

(Eburst OC R 2 ). 

I have concentrated on the recurrence time and total burst 
energy, and have not discussed other properties, such as the 
peak luminosity or burst duration. Type I bursts from 4U 1820- 
30 generally show radius expansion, in which the luminos- 
ity reaches the Eddington limit. Kuulkers et al. (2003) re- 
cently surveyed the peak luminosity of globular cluster bursters. 
For 4U 1820-30, the average peak flux is 5.27 ± 0.72 x 
10" 8 erg cm" 2 s" 1 , or a luminosity 3.7 x 10 38 erg s" 1 for a dis- 
tance of 7.6 kpc. The pure helium Eddington limit is 3.5- 
5 x 10 38 erg s" 1 for M = 1 .4-2 M Q (not including general rela- 
tivistic corrections). When hydrogen is present, the Eddington 
luminosity is reduced by a factor 1 + X. However, for small 
X, this change is much less than the uncertainty in the neutron 
star mass, and so does not allow a determination of hydrogen 
content. 

The durations of bursts from4U 1820-30 are « 20-30 s (with 
exponential decay times < 10 s), consistent with pure helium 
burning, in which the helium burns rapidly and the burst du- 
ration is set by the cooling time of the layer. For solar abun- 
dance of hydrogen, slow hydrogen burning via the rp process 
(Wallace & Woosley 1981) extends the decay time to ~ 100 s 
(Hanawa, Sugimoto, & Hashimoto 1983; Wallace & Woosley 
1984; Hanawa & Fujimoto 1984; Schatz et al. 1998). However, 
forXo < 0.1, most of the protons are able to burn by direct cap- 
ture on carbon. Therefore the small amount of hydrogen con- 
sidered in this paper is not expected to substantially increase 
the burst duration. This is something that could be addressed 
by time-dependent simulations with increasing hydrogen frac- 
tions. 

3.2. Superbursts 

In §3, I used the ra's determined by matching the Type I 
burst recurrence times to calculate superburst ignition condi- 
tions. For pure He accretion, the expected superburst recur- 
rence times are w 1-2 years (for (2 crust ~ 0.1-0.2 MeV per nu- 
cleon), much less than found by SB02 (who took m = m x ). If 
hydrogen is present, the smaller m gives recurrence times « 5- 
10 years. This approach of simultaneously modelling Type I 
bursts and superbursts can be usefully applied to other super- 
burst sources which exhibit regular Type I bursts. A promising 
candidate is KS 1731-260, for which observations of the qui- 
escent flux constrain Q cmst (Rutledge et al. 2002). In addition, 
these results emphasize that determining the superburst recur- 
rence time would strongly constrain the local accretion rate and 
thermal structure of the star, and thereby composition of the 
accreted material. 
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A direct effect of hydrogen in the accreted material will be 
to reduce the amount of carbon made during H/He burning. 
For pure helium accretion, SB02 argued that carbon produc- 
tion occurs during stable burning at high m, for which the tem- 
perature is low enough to give a high carbon yield (Brown & 
Bildsten 1998). However, if hydrogen is present, protons will 
rapidly capture on 12 C, initiating the hot CNO cycle. In that 
case, breakout reactions 15 0(a,7) and 14 0(a,p) become possi- 
ble (Wallace & Woosley 1981; Schatz et al. 1998), depleting the 
CNO abundance, and reducing the carbon yield. For Xq = 0.1, 
there are enough protons for one proton capture on each 12 C. 
The outcome is probably different for X c above and below this 
value, which is within the range predicted by models (Podsiad- 
lowski et al. 2002). For solar abundance of hydrogen, Schatz 
et al. (1999) found stable burning gave Xc ~ 10% (still enough 
to ignite a superburst, however, as CB01 showed). Calculations 
of the carbon production for stable helium burning with a small 
amount of hydrogen should be carried out. 

Finally, it has been noted that the superburst from 4U 1820- 
30 is different from the superburst seen in other sources (Ku- 
ulkers et al. 2002). It was the most energetic, and is the only 
superburst so far to reach Eddington luminosity. This is proba- 
bly due to the low hydrogen abundance in the accreted material. 
A lower hydrogen abundance allows greater carbon production, 
and makes lighter ashes (the heavy nuclei are probably iron 
group with A « 56 rather than rp process ashes with A ~ 100) 
reducing the opacity and giving a deeper ignition. Both these 
effects give a larger nuclear energy release, which could explain 
the more energetic and luminous superburst. 

3.3. Transition to Stable Burning 

The calculations described here apply only to the low state, 
when regular Type I bursts are seen, and do not address the 
disappearance of bursts in the high state. This is an impor- 
tant question to answer. In spherically-symmetric models, 
burning becomes stable at high m's because the temperature- 
sensitivity of the triple alpha reaction decreases with tempera- 
ture, and the incoming helium then burns at the rate at which 
it accretes (Fujimoto, Hanawa, & Miyaji 1981; Ayasli & Joss 
1982; Taam, Woosley, & Lamb 1996). However, the tran- 
sition to stable burning for pure helium accretion occurs at 
m st ab ~ 2 x 10 6 g cm -2 s _1 (B95; Bildsten 1998), much higher 
than 1820-30's accretion rate. The extra nuclear energy release 
if hydrogen is present will decrease m sta b, but not enough to 



explain the observed transition. B95 suggested that the miss- 
ing piece of physics was that the bursts become radiative rather 
than convective in the high state. The burning then proceeds as 
a slowly propagating "fire" over the stellar surface rather than 
dramatic Type I bursts. However, the ignition of a local spot on 
the surface of the star is not understood, requiring stabilization 
of pressure gradients, perhaps by rapid rotation (Spitkovsky, 
Levin, & Ushomirsky 2002; Zingale et al. 2002). 

4U 1820-30 fits into the pattern of bursting seen in other 
sources: a transition from frequent, regular bursting to in- 
frequent, irregular bursting as accretion rate increases (van 
Paradijs, Penninx, & Lewin 1988). By comparing Bep- 
poSAX/WFC data for 9 frequent bursters (including 4U 1820- 
30), Cornelisse et al. (2003) concluded that this transition oc- 
curs at a universal luminosity Lx ~ 2 x 10 37 erg s" 1 . If this 
behavior is common to all bursters, two explanations that have 
been put forward for hydrogen accretors may be relevant for 
4U 1820-30. First, Narayan & Heyl (2003) recently studied the 
linear stability of quasi-steady burning shells, finding that some 
stable burning occurs during accumulation for accretion rates 
near the transition to stability. Whether this result applies to 
pure helium accretion should be investigated further. Second, 
for hydrogen accretors, Bildsten (2000) suggested that the frac- 
tion of the star covered with fresh fuel increases with increasing 
global M, so that the local accretion rate m decreases. The ob- 
served transition is then from mixed H/He ignition (f reC ur < fe) 
giving frequent, regular bursts to pure He ignition, giving irreg- 
ular, less frequent bursts (f reC ui > tn)- The different composition 
in 4U 1820-30 makes it difficult to apply the same explanation, 
and the regularity of the bursts perhaps implies complete cov- 
ering. Increasing area with M would explain the small change 
in recurrence time seen by Clark et al. (1977). Accurate flu- 
ence measurements together with spectral fits to the radius as a 
function of f recur are a way to test this picture. 

I thank Erik Kuulkers, Phillip Podsiadlowski, Hendrik 
Schatz and the referee for useful comments, and Lars Bildsten 
for stressing that Type I bursts might test evolutionary mod- 
els for LMXBs. This work was supported by NASA through 
Hubble Fellowship grant HF-01 138 awarded by the Space Tele- 
scope Science Institute, which is operated by the Association of 
Universities for Research in Astronomy, Inc., for NASA, under 
contract NAS 5-26555. 
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FIG. 1 . — Temperature profiles at ignition for accretion of pure He and H/He with X = 0. 1 at a rate m = 1 .2 X 10 4 g cirT 2 s _1 . The extra heating from the hot CNO 
cycle results in a significantly hotter layer, and earlier ignition. The dashed curve is the ignition curve for pure He. The model with hydrogen ignites a little below 
the ignition curve due to the enhancement of e^ a by proton captures (see text). 
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FIG. 2. — Burst energy and recurrence time for different choices of instantaneous accretion rate m. I show hydrogen models with X = 0.1 (long-dashed line) and 
X = 0.2 (solid line), and helium models with Qcmst = 0.1 MeV per nucleon (dotted line) and Qcnist = 0.2 MeV per nucleon (short-dashed line). I take (m) = 2m at 
each m, and Z = 0.01. The vertical bars show the range of observed properties, and have been placed at the m inferred from the accretion luminosity rhx- 
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FIG. 3. — Variation of burst properties with persistent flux for a pure He model (dotted line), a low metallicity hydrogen model (short-dashed line, X = 0.01, 
Z = 10~ 3 ), and two high metallicity hydrogen models (Z = 0.01; X = 0.1 long-dashed line; X = 0.2 solid line). All models have Qcrust = 0.1 MeV per nucleon. 
The fluence and recurrence time have been normalized to the values observed by Haberl et al. 1987, fluence 3.5 X 10~ 7 erg cm~ 2 and « recll r = 3.2 hours at Fx = 
4.3 X 10 -9 erg cirT 2 s -1 . The kinks in the solid and long-dashed lines occur when the recurrence time equals the time to burn all the hydrogen (when Z = 0.01, this 
is m 3 hours for X = 0.1, « 6 hours for X = 0.2). 
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FIG. 4. — Predicted burst fluence against recurrence time for the models of Figure 3. Little variation of fluence with < recll r is seen for pure He accretion (dotted 
line), or when the metallicity in the layer is low (Z = 1CT 3 , dashed line). The hydrogen models with Z = 0.01 show 10% variation in fluence across the observed 
range of recurrence times. For X = 0. 1, the slope changes at < rccU r = 3 hours, when the time to burn all the hydrogen equals the recurrence time. 



Table 1 

Ignition conditions for bursts with f recur = 3.2 hours" 



Model 


x 


z 


6crust 


rri4 




r 8 


x b 


<*> 


2nuc 


£39 C 


1 


0.0 


0.01 


0.1 


2.3 


2.7 


1.7 


0.00 


0.00 


1.6 


3.9 


2 


0.0 


0.01 


0.2 


1.7 


1.9 


1.8 


0.00 


0.00 


1.6 


2.8 


3 


0.1 


0.01 


0.1 


1.4 


1.6 


1.8 


0.00 


0.05 


1.8 


2.6 


4 


0.2 


0.01 


0.1 


1.4 


1.6 


1.9 


0.09 


0.14 


2.2 


3.2 


5 


0.3 


0.01 


0.1 


1.4 


1.6 


1.9 


0.19 


0.24 


2.6 


3.9 


6 


0.1 


0.01 


0.2 


1.2 


1.4 


1.9 


0.00 


0.05 


1.8 


2.3 


7 


0.1 


0.001 


0.1 


2.0 


2.3 


1.7 


0.09 


0.09 


2.0 


4.2 



a columns are: accreted hydrogen fraction Xq; CNO mass fraction Z; local accre- 
tion rate m 4 = m/10 4 g cm -2 s" 1 ; ignition column depth y 8 = y/10 8 g cm" 2 ; ignition 
temperature T% = T/10 8 K; hydrogen mass fraction at the base X; mean hydrogen 
mass fraction {X) = J dyX(y)/y; nuclear energy release, Q nuc = 1.6 + 4.0 {X} MeV 
per nucleon; predicted burst energy £39 = £burst/10 39 erg. 



Table 2 

Expected properties of Type I X-ray bursts from evolutionary models : 





X() m frecur (h) 


x h 


Qnuc 


E39 




Podsiadlow 


ski et al. (2000) 






8.8 


0.00 2.2 


3.5 


0.00 


1.6 


4.1 


4.3 


0.35 1.1 


4.3 


0.20 


2.7 


4.3 


11 


0.18 2.8 


1.4 


0.13 


2.2 


2.9 




Fedorova & 


:Ergma(1989) 






15 


0.00 3.6 


1.4 


0.00 


1.6 


2.6 


4.1 


0.13 1.0 


4.9 


0.00 


1.8 


3.0 


5.8 


0.11 1.5 


2.9 


0.01 


1.8 


2.6 


7.0 


0.09 1.8 


2.3 


0.01 


1.8 


2.4 


9.6 


0.05 2.4 


1.6 


0.00 


1.7 


2.2 


a given 


the time-avera 


ged accretion 


rate 


(MU = 



(M)/10 9 M Q yr 1 and hydrogen fraction at the surface 
of the secondary Xq; I assume (m) = 2m, giving rh\ = 

(M)_ 9 /4.0 for R= 10 km. 
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Table 3 

Superburst Ignition Models" 



Model b 




(m) 4 


Gcrust 


frecur (jT) 


yi2 




£43 


-Xc,min° 


SB02 




2.4 


0.1 


14 


11 


3.8 


3.9 


0.13 


1 


0.0 


4.6 


0.1 


1.9 


2.8 


5.0 


1.0 


0.08 


2 


0.0 


3.4 


0.2 


1.1 


1.2 


5.6 


0.43 


0.16 


3 


0.1 


2.8 


0.1 


9.2 


8.1 


4.0 


3.0 


0.12 


6 


0.1 


2.4 


0.2 


3.7 


2.8 


4.9 


1.0 


0.19 



a (m) 4 = (m}/10 4 grants- 1 ; yn = y ign /10 12 gem" 2 ; £43 = 
£superburst/10 43 erg. We assume a composition of 30% 12 C, 70% 56 Fe, 
and calculate the total energy release ^superburst = 47r/? 2 yi gn gnuc assuming 
2nuc = 1 MeV per nucleon. 

b The model number refers to Table 1, except for SB02 which refers to 
Strohmayer & Brown 2002. Models 1 and 2 are pure He accretion; models 
3 and 6 have X = 0.1. 

The minimum mass fraction of carbon required for unstable ignition. The 
carbon burns stably during accumulation for lower X c . 



